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Abstract 

Lunar meteoroid impact flashes provide a method to estimate the flux of the large 
meteoroid flux and thus their hazard to spacecraft. Although meteoroid impacts on the Moon 
have been detected using video methods for over a decade, the difficulty of manually searching 
hours of video for the rare, extremely brief impact flashes has discouraged the technique's 
systematic implementation. A prototype has been developed for the purpose of automatically 
searching Lunar video records for impact flashes, eliminating false detections, editing the 
returned possible flashes, and archiving and documenting the results. The theory and 
organization of the program is discussed with emphasis on the filtering out of several classes of 
false detections and retaining the brief portions of the raw video necessary for in depth analysis 
of the flashes detected. Several utilities for measurement, analysis, and location of the flashes 
on the moon included in the program are demonstrated. Application of the program to a year's 
worth of Lunar observations is discussed along with examples of impact flashes as well as 
several classes of false impact flashes. 
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Lunar Flash Video Search - Method 
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• Detection based on assumed 
instrument point spread function (PSF) 

• Complete PSF pattern must exceed 
chosen local statistics 

• At least one pixel in the PSF must 
exceed a chosen peak sigma 

• Lunar masking & dark sky masking 
used to limit search and define 
statistical parameters 

• Prototype software, LunaCon, was 
developed by W. Swift in IDL 

• Current software, LunarScan, 
independently developed by Peter 
Gural, is fast and flexible. 


Pixilated point spread function search patterns 
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Lunar Flash Video Search - Software 


Full Frame 


Frame : 



Flash Editor (LunaCon) 
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^Software is semi-automated. Possible flashes need to be edited by operator 
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Lunar Flash Video Search - Artifacts 


I tinar f Unit Suspect ( dHor G U unar Fideo\May?006\May7*vi 


€J I unar FUth Suspect I ditor F:« unar V<droUuly?006Uuty? av» 


| Frame : 

96647 

Center : 

(565. 174) 

1 [Field : 

193295 

Peak Sigma: 7.48 

1 |lnten: 

147.19 

LAi : 

1.00 

: ««i 

< | Use 

J > 

>» | Close | 


Flash Index, Nflash - 9 



Frame: 


6753 Center: (475. 389) 


±1 


[Field: 

13507 

Peak Sigma: 2S.72 

1 Inten: 

749.47 

L/W : 1.00 

«< | 

< | Use 

>|>>>| Close | 

13 


Flash Index. Nflash - 31 


Frame 96653.0 Inten - 94.0000 Ulh - 19.4458 | 19:26:44 

. . * * ■ T’ -r/". r. • " ». 

j*£~* ■ 'j. L - 


3 Lunar Flash Suspect Editor G:\Lunar_VideoVMay2006VMay2.avi 



Fromes after 964*4 


Noise 


ii 


♦000 E 


>9 1 unar Flash Suspect f ditor G U unar Video\June7006\4June.avi 



|f~F 

118227 Center: (105, 143) 

| [Field 

236454 Peak Sigma: 44.89 

r" 

1735.15 L>W : 1.00 

< 1 Use If - > | »> | Clos 


Jj 


22 

J 

Flash Index. Nflash - 74 


*. Ji fc fr %\ t* b L 4- 



4 6 

Framat after 116224 




■Frame 

[Field 


I |lnten 


92933 Center: (609, 326) 


185867 Peak Sigma: 76.46 


<<< 


1496.16 L^W : 2.94 

< I Use !>!>>> 


Close 


8 


Flash Index, Nflash - 15 




Interpolated Field Light Curve 


Frames after 6750 
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Lunar Flash Video Measurement - Lunar Area 


• The Earth Shine area of the moon varies with phase and position. 

• Area observed must be accumulated, pixel by pixel, throughout the observations. Lunar disc 
center pixel area, A c , km 2 is determined by Lunar disc radius in pixels. 

• Surveyed area is weighted by pixel distance, r/R, from center of the Lunar disc to account for 
spherical Moon effects. 

• Lunar disc center pixel area, A c , km 2 is determined by Lunar disc radius in pixels. 

• Observed Area-Time product (km 2 -hr) ,AT, is binned by “limiting magnitude”. 

• Good magnitude adjusts for clouds, dropouts, and the “limiting magnitude”. 

• Area-time product of null observations is part of flux estimate. 

• Example: 355mm diameter f/8 telescope with 0.33 focal reducer operating at f = 950mm 

Surveys ~50km 2 -hr / frame or 1 hour video => 5e6 km 2 -hr 
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Lunar Flash Video Measurement - Photometry 


Problem 

• One has no choice of sky conditions or calibration stars 

• Available star can be hours away and at different elevation 

• Sky conditions can vary from frame to frame 

• Observations are often continued to very low elevations with 
rapidly changing airmass 

Solution 

• Use mean lunar intensity, Lmi, as transfer standard 

• Lunar sensitivity function, Lsf, includes extinction 

• SenO k can be determined for every field k , IU 

• Flash light curve, i k , IU, using Circular Aperture Photometry 

• Calibration using SenO k removes atmospheric effects 


Lsf = 


Sen 0 


o 


Apsf 0 ( Lmi 0 - Smi 0 ) 


SenO k = Lsf x Apsf 0 x [Lmi k - Smi k ) 
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Lunar Flash Video Measurement- Moon 



Definitions: 

IU 

Instrument Units 

Lmi 

Lunar mean intensity, IU 

Smi 

Sky mean intensity, IU 

[Max, Min] 

95 th and 5 th percentiles 

PSF 

Point Spread Function sigma 

Apsf 

Area of PSF, 7 r(PSF) 2 , pixels 

Lpsfi 

Lunar PSF intensity, IU 

Lpsfm 

Lunar PSF intensity, mag. 

SenO 

Sensitivity, senO = IU m 0 star 

C 

Contrast of lunar image [0,1] 

range 

mean ( 1± C), see [Max, Min] 

Range Lpsfm 

Lunar PSF magnitude range 

Lsf 

Lunar sensitivity factor 


^ _ Max - Min 

Contrast , C = 

Max + Min 


range = Lmi(1±C) 


Lpsfm image =-2.5Log XQ 
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\~5Log XQ 

PSF ,„r j 


Range Lpsfh, 


-2.5 Log xo 


' Lpsfi 
^ senO j 


+ 


[-2.5Log l0 (l + C),-2.5Log l0 (\-C)] 


>Range Lpsfm can be used as proxy for Limiting Magnitude for Highlands and Mare 
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Lunar Flash Video Measurement 
Flash Characterization 


• l Totaj and maximum intensity, l a often the only measurables of Lunar flash 

• Measures are required that are independent of camera, exposure and sky 

• Most flashes can be well approximated by exponential thermal decay curve 

• Thermal Decay is specified by Initial intensity, i 0 and Time Constant, a 

• For single frame flashes, l Total = l a = i 0 and a = 1 (by definition) 

• Flash Initial magnitude m f and Time Constant a form ideal measures 

• m f — 0 — > i 0 = 5x1 0‘ 9 watts/m 2 so: m f = 6, a = 0.10 sec => 2x10 -12 J/m 2 



Time Constant, a = 


a 




Initial Intensity, 



Flash Initial magnitude, m f = -2.5 Log l0 




\ 




sen 0 


k J 
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Lunar Flash Video Measurement 
Characterization Example 


May 2, 2006 Impact: 

Intensity Plot 

• Time Constant, a = 0.04 

• Initial Magnitude, wy= 6. 1 

• Maximum recorded, m a = 6.8 




May 2, 2006 Impact: 

Magnitude Plot 

• Decay linear in magnitude 

• Slope yields the time constant 

• Very bright flash can exhibit 

complex behavior. 


Seconds 



Data 


Analysis 
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Lunar Flash Archiving - 
Reporting 



www.nasa.gov/centers/marshall/news/lunar/index.html 
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Summary 




• Methods have been developed for semi-automated lunar impact flash detection. 

• Software and techniques have been developed to edit suspected lunar impact flashes. 

• Methods are described for finding Lunar survey area as a function of detection limits. 

• Lunar intensity over the PSF is developed as a photometric calibration transfer standard. 

• A system for characterization of impact flashes is developed independent of instrumentation. 

• Lunar flash data archiving techniques are described and preliminary survey results shown. 
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